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C, of the tangent, at the point P, to the given helix. This tangent cuts the tan- 
gent t', which is drawn at the point P' to the circle T, in the horizontal trace T\ 
of the required tangent. Etc. 
Genoa, May-June, 1918. 

DIFFERENCE QUOTIENTS. 

By J. P. BALLANTINE, Cambridge, Mass. 

The whole subject of interpolation, in the case of equally spaced ordinates, 
by means of algebraic polynomials is commonly based on the subject of finite 
differences. The method, however, fails for unequal spacing. The purpose of 
this paper is to show, that by a suitable extension of the conception of the dif- 
ference quotient, the method can be generalized for the case when the ordinates 
are spaced according to any law. 1 

Let f(x) be any single valued function. Plot the curve y = f(x). Pick out 
on the curve the (n + 1) points whose x coordinates are 1, 2, •••, (n + 1), 
respectively. By means of the method of finite differences, find the equation of 
the curve of the form y = P n (x) which passes through the (n + 1) points. 
For convenience, call this curve a secant of degree n, noticing that a secant of the 
first degree is an ordinary secant. It is easily verified that the nth derivative 
of this secant is precisely equal to the difference of order n of the ordinates of 
the (n + 1) points. 

Now let us take our (n + 1) points on the curve at random. We have to 
define the difference quotient of order n at these (n + 1) points in such a way 
that when the values of x u x 2 , • • •, x„+i are 1, 2, • • •, n + 1 respectively, it will 
reduce to a difference of order n. 

We can define the secant through these (n + 1) points as the curve of the 
form y = P n (x) which passes through them. By way of a more explicit notation, 
y = P n (x; 1, n + 1) will be the equation of the nth degree secant at the points 
(x h yi), (x 2> 2/2) •• • (Xn+u yn+i). It follows that 

P n (xi; 1, n + 1) = yi, l<j<n+l. 

Definition. The difference quotient of order n, taken at the (n + 1) 
points (xi, 2/1), (x 2 , y 2 ) ■■ • (a;«+i, y„+ 1) is the nth derivative of the nth degree 
secant through those points. It is denoted thus Qi, n +i: 

Q?, n+ i = ^Pn(x; l,n+D. 

Consider the special polynomial and associated equation 

Pr^i(x; 1, n + 2) - P n (x; 1, n + 1), 

<1) P„+i(a;»; 1, n + 2) - P n (xi; 1, n + 1) = y t - yi = 0, 1 < i < n + 1. 
1 Cf . A. A. Markoff, Differenzenrechnung, 1896, p. 10. 
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From the factor theorem for polynomials, we may say that: 

(2) Pn+i(x; 1, n + 2) — P n (x; 1, n + 1) = K(x — xjix — x 2 ) • • • (x — a^n). 
Now, determine K by (n + 1) differentiations: 

(3) QW +2 =K(n+l)L 
Using this value of K in (2), and transposing, we have 



(I) P n+ i(x; 1, n + 2) = P„(a;; 1, n + 1) + T^rrr; (a; - Xi){x—x % ) • • • (x—x^i). 



M 1 

(n+l)V 



In this way a secant of degree (n + 1) is expressed in terms of one of degree 
n and a difference quotient. We know how to write the equation of a secant 
of degree one or zero, and so our problem of writing the equation of a secant of 
the nth degree at (n + 1) arbitrary points is solved, except for determining the 
difference quotient. 

In the same way as we inferred the truth of (2), we may show that 

(2') P n+ i(x; 1, n + 2) - P n (x; 2, n + 2) = K{x - x 2 )(x -x 3 )---(x- Xn+ 2 ). 

It is not obvious that the two arbitrary constants must be the same, until 
one notices that in each case K must be equal to the coefficient of x n+1 in the 
polynomial, Pn+iix; 1, n + 2). 

Subtracting (2') from (2), we have 

(4) /Kixn+i— xi)(x— x 2 )(x— x 3 )---(x — x n ± 1 ) = P n (x; 2, n + 2) — P n (x; l,w+ 1). 
Differentiating n times, we obtain 

(5) K(x n+i - xj (nl) = Ql „+ 2 - Ql „ +1 . 
Then by means of the value of Qj^n+2 as obtained in (3), we have 

/tts rm+l _ Q",n+2 — Qi,n+1 

V 11 ^ Vl,n+2 — 

Xn+2 — X\ 

n+1 

We are now in a position to form a table of difference quotients. Formula 
(II) tells us, that besides forming a difference, we must divide by a denominator, 
which is simply the average of the (n + 1) intervals of x between the (n + 2) 
points at which Qi^n+2, the entry we are at the time computing, is taken. It 
may be denoted thus: 

(6) A k Xi = £ . 

Applying the definition of the difference quotient of order zero, we notice 
that the secant of order zero through the point (xtyi) is given by the equation: 
y = Vi- Without differentiation, we find : Q ° = yi . 



1919.] DIFFERENCE QUOTIENTS. 55 

The table of difference quotients may be arranged as follows: 



. A 2 a?i x 2 

A ' Xl A 2 x 2 A f* x B 

AlX3 x< 



Vl ~ ^' AO° } 

-,, — n° A "i. 2 "i> 2 Ar) 1 n 2 

y2 ~ ^ 2 ao° o l yi - 3 ^ 1>3 A0 2 o 3 



The elements x% and ?/i are assigned. 

All the elements to the left of these can be computed by relation (6) above, 
and those on the right by the following two relations: 

(7) AQt^ + ic + i = Qi+i, i+k+i Qi,t+k> 

/ox rt *+i _ Aft, j+k +i 

(8) %i+ * _ A^Xi • 

One way of writing the equation of the secant through the four points 
(zi2/i), (x 2 y 2 )-- • would be: 

(9) y = yi+-jy(x— xi) + -~ (x - xi)(x - x 2 ) + -^(2 — xi)(x-x 2 )(x-x z ), 

obtained by a repeated application of equation (I). 

In the general case the formation of a table of difference quotients involves a 
few laborious divisions; but these can be avoided by a proper choice of the 
spacing of the intervals of x. If, for instance, all of the intervals of x are unity, 
then AkXi, the average of any k of them, will also be unity, and the rule for the 
formation of the difference quotients reduces to that of the formation of a table 
of differences. In that case the formulae of Newton and Gauss are special cases 
of form (I). 

These are by no means the only important special cases. Let us first consider 
the limiting case when one point, {xj, yj), approaches the point, (a;,-, yi). The 
limiting value of Q\j would be /'(xi), if such a derivative existed, where /(a;) = y 
is the path along which the first point approached the second. The derivative 
may be regarded as a special case of the difference quotient, by extending the 
conception of the secant, as defined earlier, to take into account multiplicity. 
Whether regarded as a special case or as a limiting case, it may be used in 
formula (I) and on the right hand side of (II). Equation (II) does not hold in the 
special case when Q1^+ 2 = / (n+1) (^i)> for if all the in + 2) points are regarded 
as coincident, then the denominator is zero. In that case, however, / (n+1) (xi) 
is part of the data and does not have to be computed. 

Let us illustrate the case of two coincident points by a simple example. To 
form a polynomial which, when x has the value 1, has the value 2 and its deriva- 
tive the value 3; and when x has the value 4, has the value 5. The table of 
difference quotients is below. 
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A 2 x 


AiCE X 


2/ 


Ay Q 1 


AQ 1 Q 2 


1 


6 4 


2 

2 
5 


3 
3 1 


-2 -| 



Then 



2/ = 2 + rx - 1)3 + (a - 1) : 



(-9(i) 



In the special case when all the points are coincident, the difference quotients 
are derivatives, and the polynomial one obtains from equation (I), is exactly 
the polynomial one would obtain using the same number of terms in the Taylor 
expansion. 

Suppose we had an indefinite integral tabulated to a large number of decimal 
places, and wished to interpolate. For the sake of argument, suppose the table 
of differences would not converge rapidly enough. It often happens that we 
are able to compute the derivative of a function represented by an indefinite 
integral quite readily. It appears intuitively, that if we could in some way use 
to advantage the values of the derivative the process would converge more 
rapidly. 

In the case under supposition, we know the values of a function and its 
derivative for equally spaced values of the argument. It is as if we knew the 
value of the function for equally spaced pairs of points, the two points of each 
pair being coincident. The difference quotient at any pair of coincident points 
is given to us as the derivative, and all the other difference quotients may be 
found by formula (II). There are so many simplifications in the process, that 
I will submit an example. 

To find the value of logio 31.029971 to ten places. First we notice that it is 
easy to compute the derivative : 



d 



1 



dx logio x = - logio e = 



0.4342944819 



We tabulate a few nearby logarithms, and the corresponding values of 



0.4342944819 

, as follows: 

X 








X 


log X 


derivative 




29 
30 
31 
32 
33 


1. 46239 79979 
1. 47712 12547 
1. 49136 16938 
1. 50514 99783 
1. 51851 39399 


0. 01497 56718 
0. 01447 64827 
0. 01400 94994 
0. 01357 17026 
0. 01316 04389. 




The table may be arranged as follows: 






x y = log x 
29. 1. 46239 79979 

29. 1. 46239 79979 

30. 1. 47712 12547 


Q 1 

0. 01497 56718 
1472 32568 
1447 64827 
1424 04391 


iQ 2 

-0. 00025 24150 

24 67741 

23 60436 

23 09397 


0. 00000 56409 

1 07305 

51039 

97248 
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X 

31. 
32. 
33. 



y - log x 



Q 1 



1. 49136 16938 0. 01400 94994 
1378 82845 

1. 50514 99783 1357 17026 
1336 39616 



1. 51851 39399 



1316 04389 



*Q 2 

-0. 00022 12149 
21 65819 
20 77410 
20 35227 



*AQ 2 

0. 00001 46330 
88409 
42183 



00001 12818 
1 07305 
1 02078 
97248 
92660 
88409 
84366 



IQ 4 

-0. 00000 05513 
5227 
4830 
4588 
4251 
4043 



IAQ 1 

0. 00000 00286 
396 
243 
336 
207 



0. 00000 00858 
792 
729 
672 
621 



The entries in the column headed y are the data as copied from some table. 
Notice that the second entry is the same as the first, for they are the values of 
the function for two coincident values of x. This repetition is superfluous, 
and may be omitted. The column headed Q 1 contains the difference quotients 
at every two consecutive points. The first two values of x are coincident, so 
the first difference quotient is copied from the table of derivatives. The second 
and third values of x are a unit apart, so the difference quotient is the difference. 
And so on down the column, alternate entries are copied from the table of deriva- 
tives, and the others are obtained as differences of y. 

In order to obtain Q 2 , formula (II) tells us to divide the difference of two con- 
secutive entries in the column of first order difference quotients by J. If you 
have a set of intervals, alternate ones zero and one, the average of any two 
consecutive ones of them will be f . Let us omit dividing by this factor, and 
thereby obtain ^Q 2 instead of Q 2 . To obtain the entries in the column headed 
|AQ 2 , we difference those of the column before. Formula (II) says that in order to 
obtain Q a we must divide the first AQ 2 by 1/3, and the second by 2/3. The dif- 
ferences of x are 0, 1, 0, 1, etc., and the average of three consecutive intervals 
will be 1/3 and 2/3 alternately. In order to avoid labor, we multiply every other 
one by 2, for they look smaller than the others, and thereby in each case obtain 
|Q 3 instead of ^Q 3 . Absorb these factors into the difference quotients, denoting 
the pseudo difference quotients by Qi n . The factorial divisors in the terms of the 
polynomial will be changed, as follows: 



(10) 



**(*) = Vi + 

+ 
+ 



tQi 



+ 



PQi 2 



*»(«.- l)Qi» , fit - l) 2 Qi 4 



(0!) 2 (1!K(0!) 2 (1!) 2 ^ (1!) 2 (2!) 

f(t - l) 2 (t - 2) 2 - • • ft - n + l) 2 (t - n) 
(n\) 2 (n+ 1)! 

t 2 (t- l)\t-2) 2 ---{t-n) 2 



(1!)(2!) 2 



+ 



Qi 



2»+l 



(n!)((n+l)!) 2 



Qi M , 



(t = x — x{). 
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We obtain in this particular example the polynomial : 
y = 1.4913616938 + 0.0140094994* - 0.000221212149* 2 

y2/y _ 1\ j2(4 _ 1\2 

+ -~ — - 0.0000092660 - 0.0000004251 -\ . 

To obtain log 31.029971, set t = 0.029971; adding the various terms of the 
polynomial, we obtain the result 

1. 49136 16938 

41 98787 

-1995 

—41 

-1 



1. 49178 13688= log 31.029971. 

Although the definition given above for the secant is geometrical, it could 
just as well be stated analytically. The work which I have carried out for reals 
could be extended into the domain of imaginaries. 

When one has obtained the interpolatory polynomial, the secant, it may 
be asked, "How near will this come to giving the true value of the function 
throughout a certain interval?" 

In equation (I), set x = Xn+2, and replace Pn+z(xn^2; 1, n + 2) by its 
equivalent value /(a; n+2 ). Considering a;„+ 2 as the independent variable, and 
calling it x for simplicity, we get the equation: 

(III) f(x) =- P n (x; 1, n + 1) + { Q?X\i(x) } (* - *i) (x-x 2 )---(x- x^) ■ l 



(n+l)l 

where Q™1»+i(a;) is the difference quotient taken for the (n + 1) fixed values 
#i, #2, • • • , Xn+i, and the one variable value, x. It is a function of x. 

A remainder is that which has to be added to an approximation in order to 
give the true result. 

Let 

f(x) = P n (x; 1, n + 1) + R(x). 
Then 

(IV) R(X)= , { Q?+Vl (X) }(X - X!)(X - X 2 )- ■ -(X - Xn+l). 

We have determined R(x) in terms of Qffl+iix), and so let us study the latter. 
Let us examine its derivative. By the use of formula (II), we have 

<*Vl,n+lW J __ -r . Ql,n+2 Qo,n+l _ j . Qo, n+2 



UX \x=x aJ„ + 2-»«o x n+2 Xq a:„«-»a-o n "t" " 

In general there is no limit, but let us suppose f'(x) exists at x = x - Then 
the secant at the points Xo, X\, • • • , x„ +2 approaches a limiting secant as a;^.2 
approaches x , namely the secant at the points x , x\, • • • , Xn+i, which has contact 
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at Xo. Its derivative of order (n + 2), the difference quotient, has a limiting 
value also, and is denoted thus: Qo,"o?»+2- Then we have 



(V) 



dQVi+t(x) 



dx 



n + 2 



Qn+2 



n+2 



n + 2 



Proceeding in this way one can determine the properties of Q"+J +1 (x) from 
those of fix), and in that way can obtain the oscillation of the remainder for 
any interval. If the problem is to determine the derivative of fix), one differ- 
entiates in formula (III), and obtains the expression for the remainder in that 
case in terms of Qffl+i(x) and its derivative, which are already known. 

In this way, the difference quotient enables us to expand an arbitrary function, 
analytic or merely continuous, real or complex, at points to suit the peculiar 
needs of the problem. Moreover, in any special case, in which a simpler and less 
elastic method of expansion can be used, the expansion in terms of difference 
quotients reduces automatically to that method. The remainder can be ex- 
pressed explicitly in terms of the difference quotient, not as in the case of a 
remainder in terms of a derivative at some unknown point, but taken at the 
points about which we are expanding and the one additional point x, the point 
at which we want the remainder. 



A THEOREM IN THE GEOMETRY OF THE TRIANGLE. 

By J. W. CLAWSON, Ursinus College. 

It is the purpose of this paper to state a general theorem of which the well- 
known Wallace (or Simson 1 ) line theorem is a very special case; to point out the 
result of a polar reciprocation; to mention some special cases of the theorem; 
and to give an application to the geometry of the quadrilateral. 

1. Theorem? AiA^Az is a triangle, a< the side opposite Ai, (i = 1, 2, 3), 
P a point on the circumference of the circumscribed circle 3 C, Q any point in the 
plane. Let AiQ cut the circle 3 at J7,; PUi intersects a» at Li. Then L\, L 2 , L 3 , Q 
are collinear on I. Let us call this line the "P-transversal of Q" with respect to 
the triangle A1A2A3. 

The theorem follows at once from Pascal's Theorem for, since A2A3A1U1PU2, 
A1A2A3U3PU1 are inscribed hexagons, L u X 2 , Q and L z , L u Q are collinear. 

Cor. 1. — If a transversal cut the side a» of a triangle at Li, and P be any 
point on the circumcircle, and if PLi cut the circle at £/*, then A1U1, A2U2, A3U3 
are concurrent at Q, a point on the transversal. 

Cor. 2. — If a transversal cut the side a t of a triangle at Li, if Q is any point 

1 For references, see this Monthly, 1916, p. 61. 

2 It would be strange if this simple theorem were new; yet the writer has not found it 
recorded. A special case is given in Hatton's Projective Geometry, 1913, Ex. p. 164. 

3 Or conic. 



